PLANE WAVEGUIDES WITH CORNERS IN THE 
SMALL ANGLE LIMIT 

MONIQUE DAUGE AND NICOLAS RAYMOND 

ABSTRACT. The plane waveguides with corners considered here are infinite V-shaped strips 
with constant thickness. They are parametrized by their sole opening angle. We study the 
eigenpairs of the Dirichlet Laplacian in such domains when this angle tends to 0. We provide 
multi-scale asymptotics for eigenpairs associated with the lowest eigenvalues. For this, we 
investigate the eigenpairs of a one-dimensional model which can be viewed as their Born- 
Oppenheimer approximation. We also investigate the Dirichlet Laplacian on triangles with 
sharp angles. The eigenvalue asymptotics involve powers of the cube root of the angle, while 
the eigenvector asymptotics include simultaneously two scales in the triangular part, and one 
scale in the straight part of the guides. 



1. Introduction and main results 

1.1. Motivations. Quantum waveguides refer to meso- or nanoscale wires (or thin sheets) 
inside electronic devices. They can be modelled by one-electron Schrodinger operators with 
potentials having high contrast in their values. In many situations, such Schrodinger operators 
can be approximated by a simple Laplace operator with Dirichlet conditions on the boundary 
of the wires [15]. The presence of bound states is an undesirable effect which is nevertheless 
frequent and useful to predict. The same Laplace-Dirichlet problems arise for TE (transverse 
electric) modes in electromagnetic waveguides [9]. 

This is a well-known fact, from the papers [17, 15, 10, 11], that curvature makes discrete 
spectrum to appear in waveguides. Moreover the analysis of this spectrum can be accurately 
performed in the thin tube limit (in dimension 2 and 3, see [15, Section 5]). In fact, this 
asymptotical regime corresponds to a semiclassical limit so that the standard techniques of 
[24] could have been used to investigate that problem. 

Since curvature induces discrete spectrum, this is then a natural question to ask what hap- 
pens in dimension 2 when there is a corner (which corresponds to infinite curvature): Does 
discrete spectrum always exist in this case? This question is investigated with the L-shape 
waveguide in [18] where the existence of discrete spectrum is proved. For an arbitrary angle 
too, this existence is proved in [3] and an asymptotic study of the ground energy is done when 
9 goes to | (where 9 is the semi-opening of the waveguide). Another question which arises is 
the estimate of the lowest eigenvalues in the regime 9 — > 0. This problem is analyzed in [9] 
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through matched asymptotic expansions and electromagnetic experiments. This is precisely 
the question we tackle in this paper: We are going to prove rigorously complete asymptotic 
expansions for the eigenpairs in plane waveguides with corner (also called "broken strips") as 
9 tends to 0. We have provided in [14] numerical experiments by the finite element method for 
this situation too. 

For the case of dimension 3, we can cite the paper [19] which deals with the Dirichlet 
Laplacian in a conical layer. In this case, there is an infinite number of eigenvalues below the 
essential spectrum. The other initial motivation for the present investigation is our previous 
work [4] in which we study the Neumann realization on R 2 ^ = {(s, t) G IR 2 : t > 0} of the 
Schrodinger operator — d 2 — d\ + (t cos 9 — s sin 9) 2 in the regime 9 — > (see also [28, 25]). It 
turns out that the lowest eigenfunctions of this operator are concentrated near the cancellation 
line of the potential, which also enlighten the link between a confining electric potential and a 
strip with Dirichlet boundary conditions. 

In our way towards the analysis of plane waveguides with corners, a natural step turns out 
to be the study of the Dirichlet problem on isosceles triangles with small angle. This subject 
is already dealt with in [20, Theorem 1] where four-term asymptotics is proved for the first 
eigenvalue, whereas a three-term asymptotics for the second eigenvalue is provided in [20, 
Section 2]. In fact the spectral analysis of triangles with small angles is not the sole way 
to succeed in the study of waveguides. Nevertheless, as just mentioned, this problem has a 
particular interest on its own and permits to enlighten the presentation of the proofs. 

Finally, in the same vein, we can mention the papers [21, 22] whose results provide two- 
term asymptotics for the thin rhombi and also [5] which deals with a regular case (thin ellipse 
for instance), see also [6]. 

1 .2. The Dirichlet Laplacian on the broken guide. Here we introduce the family of broken 
guides tt e , parametrized by the angle 9, and give basic properties of the spectrum of the positive 
Laplacian with Dirichlet condition in £l e . Then we state our main result related to the behavior 
as 9 — > of the lowest eigenvalues of these operators. 

1.2.1. Basic properties. Let us denote by (x±, x 2 ) the Cartesian coordinates of the plane and 
by = (0, 0) the origin. The positive Laplace operator is given by — dj — d\. The domains 
of interest are the "broken waveguides" which are infinite V-shaped open sets: For any angle 
9 G (0, |) we introduce 

(1.1) Q = \(xi,x 2 ) G R 2 : Xitantf < \x 2 \ < \Xi + -7^) tan#) . 



Note that its width is independent from 9, normalized to it, and 9 represents the (half) opening 
of the V, see Fig. 1. The limit case where 9 = | corresponds to the straight strip (— n, 0) x R. 
The aim of this paper is the investigation of the lowest eigenvalues of the positive Dirichlet 
Laplacian A®" in the small angle limit 9 — > 0. 

The operator A^' r is a positive unbounded self-adjoint operator with domain 




Dom(A^ r ) = {i, G Hl$l e ) ■ ^ e L 2 (n e )}. 
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Figure 1 . The broken guide £lg (here 9 = |). Cartesian and polar coordinates. 

When 9 = f, we simply have Dom(Ag; r ) = H 2 n H%(£1 ). In contrast, when 9 G (0, f ), the 
boundary of Q e is not smooth, it is polygonal. The presence of the non-convex corner with 
vertex is the reason for the space Dom(A^ r ) to be distinct from H 2 D H^flg). Nevertheless 
this domain can be precisely characterized as follows. Let us introduce polar coordinates (p, (p) 
centered at the origin, with ip = coinciding with the upper part x 2 = %i tan 9 of the boundary 
of VLq. Let x be a smooth radial cutoff function with support in the region x\ tan 9 < |a?2 1 and 
X = 1 in a neighborhood of the origin. We introduce the explicit singular function 

(1.2) ^ ng (x ll x 2 )= X (p)p v/u ^—, with u = 2(ir-6). 

Then there holds, see the classical references [27, 23]: 

(1-3) Dom(A^) = {H 2 n © [^Jj 

where f ifing] denotes the space generated by ^fi ng - 

We denote by p>Gu\,n(9) its n-th Rayleigh quotient, n > 1 (here || ■ || is the L 2 norm on fig): 



IIVVII 
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^Gm,n{9) = inf SUp 

ipi,...,ipj independent in H^(Q e ) i>&pim{'tp 1 ,...,ipj} 

We gather in the following statement several important preliminary properties for the spectrum 
of A®". All these results are proved in the literature. We briefly indicate hereafter what are the 
main arguments of the proofs, and where details can be found. 

Proposition 1.1. (i) If 9 = |, A^ r has no discrete spectrum. Its essential spectrum is the 
closed interval [1, +oo). 

(ii) For any 9 G (0, |), the essential spectrum of A®" coincides with [1, +oo). 

(iii) For any 9 G (0, |), the discrete spectrum of A®" is nonempty and finite. In other words, 
A®" has at least one eigenvalue below 1, but a finite number of them. 

(iv) For any 9 G (0, |) and any eigenvalue in the discrete spectrum of A®", the associated 
eigenvectors if; are even with respect to the horizontal axis: ip(xi, — x 2 ) = ip(%i, x 2 ). 



(v) For any n > 1, the function 9 \-> pQ U \^ n {9) is continuous and non decreasing on (0, f ). 
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(vi) For any n > 1 and 9 such that /UG u i,n(#o) < 1» the function 9 h-> /iGui,n(0) zs strictly 
increasing on (0, # ]. 

Proof, (i) is a clear consequence of the separation of variables in VL v i 2 = (— 7r, 0) x R. 

(ii) is a consequence of the fact that outside a compact set, fig is the union of two strips 
isometric to (0, +oo) x (0, n). 

(iii) The fact that there are eigenvalues below the essential spectrum is known since [3]. See 
also in [14, §4] another proof based on a more general argument developed in [15, 10, 11] for 
waveguides with curvature. The fact that there is only a finite number of such eigenvalues is 
proved in [14, §5] using a similar method as [30, Theorem 2.1]. 

(iv) Since the domain and the operator are invariant by the symmetry x 2 — £2, the eigen- 
vectors are even of odd with respect to the horizontal axis. An argument of monotonicity for 
Dirichlet eigenvalues excludes the odd eigenvectors, see [14, §2.2] for details. 

(v) The Rayleigh quotients are non-decreasing functions of 9 as a consequence of the pre- 
vious point and a suitable change of variable which transform the operator —A in a domain 
depending on 9 into an operator depending on 9 on a fixed domain, see [14, §3] for details. 

(vi) If /UGui,n(#o) < h by points (v) and (ii), fiGu\,n(9) is an eigenvalue for all 9 e (0, 9 }. The 
same proof as in point (v) then shows that fiGu\,n(9) depend in an analytic way from 9 in (0, 9q\ . 
In addition, anticipating the result of Theorem 1.2, we find that the function 9 fJ,Q u \ in (9) is 
not constant so that we deduce from (v) that it is strictly increasing where it is analytic. □ 

1.2.2. Statement of the main result. One of the main results of this paper is a complete as- 
ymptotic expansion 1 of the eigenvalues p,cu\,n(9) in powers of 9 1 ^ 3 . To state this result, we 
need the following notation: For n > 1, let z^{n) be the n-th zero of the inverse Airy function 
A(a;) = Ai(-s). 

Theorem 1.2. For all N , there exists 9 > 0, such that for all 9 G (0, 9 ], A^ r has at least 
N eigenvalues. These eigenvalues admit the expansions: 

IM«*M ~ E%>^ /3 with = 7, 7i A n = 0, and 7 2 A n . = 2(4vrv / 2)- 2/3 ^ A (n) 

and the term of order j = 3 is not zero. The corresponding eigenvectors have multi-scale 
expansions (see Section 6.3.1 for details). 

1.3. Related questions. In the small angle limit the vertical line x± — appears as a right bar- 
rier for eigenmodes, cf. the computations in [14, §8]. In a first approach, this can be explained 
by a one-dimensional approximation in the spirit of the Born-Oppenheimer approximation: It 
is obtained by replacing — &1 2 in the expression of A^ r by its lowest eigenvalue on each slice 
of VLq at fixed x±. The effective potential thus obtained has a triangular well at X\ = (on the 
left) and a barrier on the right. That is why it is quite natural to study first a similar ID model 
operator (see Section 3). The main interest is to exhibit for such a simple situation how the 

1 By the notation X(6) ^ Sj>o c i^ P (with a positive p) we mean that for any positive integer J we have 
the estimate 

\K°) - Eo<,<,7 c j° 3P \ < C J 6> (J+1)p for 6 small enough. 
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zeros of the Airy function come into play and how two distinct scales are necessary to describe 
eigenmodes. Moreover, as a by-product of our proofs, it turns out that the first two terms in 
the eigenvalue asymptotics for A°' r and its ID approximation coincide. 




FIGURE 2. Broken guide Q g with associated triangle Trig and sector Sg. 



In the proof of Theorem 1 .2, we will have to perform an accurate analysis of the spectral gap 
separating the eigenvalues. This gap, as stated in Theorem 1.2, is of order 6> 2 / 3 and is related 
to the difference between zeros of the Airy function (in other words the gap is determined as 
soon as we have proved a two-term expansion). In order to succeed in the investigation, we 
will have to estimate this gap by comparing with a simpler spectral problem. Here we have 
to choose between several possibilities: Either we could compare with the spectrum of the 
isosceles triangle Trie (with Dirichlet conditions), or we could compare with the spectrum of 
the sector Sg (with Dirichlet conditions), see Fig. 2. The case of the sector is well-known in the 
small angle limit (we find that the first two terms in the expansion of the eigenvalues coincide 
with that of HG U \,n{Q), see [16, 20]). Nevertheless, we have preferred to analyze the problem of 
the triangle which is less known and which has an interest on its own (see [20]). In addition, 
as it will be seen in the analysis, the reduction to the triangle (through estimates of Agmon 
type) is slightly easier. A posteriori, the first two terms of the eigenvalues are the same as for 
the sector. Finally, the option to provide a full investigation of the triangle permits to divide 
difficulties inherent to each problem. This pedagogic perspective is also one of the motivations 
to study a ID model operator which roughly describes the spectral behavior of the waveguide. 

1.4. Organization of the paper. In Section 2, we discuss the different reductions to sim- 
plified operators and introduce the main notation used in this paper. We state all our results 
related to eigenvalue asymptotics. In Section 3 we investigate through a construction of quasi- 
modes and an ODE analysis the one dimensional toy model — n 2 d^ + W with the discontinuous 
triangular potential W equal to —z when z < and 1 when z > 0. In Section 4 we study a 
one dimensional approximation of the Dirichlet problem on a triangle with small angle. By 
Agmon estimates and a projection method, this leads in Section 5 to results on triangles in the 
small angle limit. Finally, in Section 6, we perform a construction of quasimodes adapted to 
waveguides and introduce in particular Dirichlet-to-Neumann operators to solve a transmis- 
sion problem; we complete the proof by comparing with the triangle case. We conclude our 
paper by discussing relations between the eigenvector asymptotics and the reentrant corner 



6 MONIQUE DAUGE AND NICOLAS RAYMOND 

singularity. We also discuss the extension of our results on X-shaped waveguides (crossing 
straight wires). 

2. Reductions 

This section is devoted to the introduction of reduced and simplified operators that we will 
consider throughout this paper. First we will use the symmetry of the waveguide to reduce 
the investigation to an half-guide. This first simplification makes a discontinuity in boundary 
conditions to appear at the origin (see Figure 3). In fact, as will be seen later, this jump in 
boundary conditions traps the eigenfunctions, which are localized in the left part of the guide. 
Due to this localization, it makes sense to tackle the Dirichlet Laplacian on triangles Trig. We 
also introduce a ID approximation of Born-Oppenheimer type for the guides and the triangles. 
This helps to understand the concentration of eigenfunctions near the origin. Finally we state 
our results concerning eigenvalue asymptotics for all these model problems. 

2.1. Half-guide and triangles. 

2.1.1. The half-guide. As a consequence of the parity properties of the eigenvectors of A^ r , 
cf. point (iv) of Proposition 1.1, we can reduce the spectral problem to the half-guide 

(2.1) Qj = {(x 1 ,x 2 ) eQ e - x 2 >0}. 

We define the Dirichlet part of the boundary by dow^g = <9f2g R <9fij~, and the corresponding 
variational space (the form domain) 

#Mix(^ + ) = t<Pe H\n+) : ^ = on 5 Dir fi+}. 

Then the new operator of interest, denoted by A^| x , is the Laplacian with mixed Dirichlet- 

Neumann conditions on . Its domain is: 

Dom(AjJ| x ) = {ipe H^ x (ttj) : Aip e L 2 (fi+) and «9 2 ^ = on x 2 = 0}. 

Then the operators A°' r and A^} x have the same eigenvalues below 1 and the eigenvectors of 
the latter are the restriction to fig of the former. 

2.1.2. Reseating of the half-guide. In order to analyze the asymptotics 9 — > 0, it is useful to 
rescale the integration domain and transfer the dependence on 9 into the coefficients of the 
operator. For this reason, let us perform the following linear change of coordinates: 

(2.2) x — X\ sin^, y = x 2 cos 9, 

which maps fig onto 0+ 4 which will serve as reference domain, see Fig. 3. That is why we 
set for simplicity 

(2.3) Q := fi+ /4 , d D -M = <9 Dir fi+ /4 , and H^ ]x (Q) = e H\n) : ^ = on <9 Dir fi}. 
Then, A^| x is unitarily equivalent to the operator defined on f2 by: 

(2.4) V Gul (9) := -2sin 2 9d 2 x -2cos 2 9d 2 y , 
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with Neumann condition on y = and Dirichlet everywhere else on the boundary of fl We 
let h = tan 9 ; after a division by 2 cos 2 9, we get the new operator: 

(2.5) C Gu] (h) = -h 2 dl-dl 
with domain: 

Dom(£ Gui (/i)) = G i^ ix (0) : C Gui {h)^ e L 2 (fi) and ^ = 0oni, = 0}. 

2.1.3. The triangles. We will also need to introduce the triangular end of this waveguide: 

(2.6) Trig = {(x 1 ,x 2 ) G Q e : xi < 0} 

and the corresponding Dirichlet Laplacian denoted by A^f . 

Prior to the investigation of £c u i(^), we are to going to study Cj r \{h) which denotes the 
same operator —h 2 d 2 — d 2 with Dirichlet conditions on the triangular end Tri of the model 
waveguide 1^/4 

(2.7) Tri = el 2 : -nV2 < x < and \y\ < x + nV2^ . 

2.2. Born-Oppenheimer approximation and models. As mentioned at the beginning of this 
section, we will use a projection method to analyze CQ U \(h). This method is based on the 
original idea of Born and Oppenheimer (see [7]) which was used to study the Hamiltonian of 
molecules (see [12, 29, 26]). By analogy with this situation, we can say that, in this paper, 
x plays the role of the nuclei variables whereas y plays the role of the electrons and where h 
would represent a mass ratio. The variable x is sometimes said to be the slow variable and y 
the fast variable. Therefore we will broaden the "molecular idea" to our waveguide situation. 

2.2.1. Schrddinger operators in one dimension. In the analysis of Cj r \(h) and £q u \(}i), we 
will see that its so-called Born-Oppenheimer approximation will play an important role: 



(2.8a) 



HmMh) = -h 2 d 2 x + V{x), 



8 



MONIQUE DAUGE AND NICOLAS RAYMOND 



where 

7T 2 

whenx G (— 7iV2,0), 



2 



(2.8b) V{x) = \ + 

- when x > 0. 

2 

Remark 2.1. This "approximation" will be justified afterwards and will not be directly used 
in the investigation. Nevertheless it already gives a non trivial insight of some fine structures 
appearing in the analysis (such as the different scalings and the transmission problem between 
the left and right parts of the waveguide). 

This effective potential V is obtained by replacing — d 2 in the expression of Cq u \{K) by its 
lowest eigenvalue on each slice of VL at fixed x. When h goes to zero, the behavior of the 
ground eigenpairs of "HBO,Gui(^) is driven by the structure of the potential near its minimum, 
attained at x = 0: In a neighborhood of x = 0, V can be approximated by its left and right 
tangents, which provides the approximate potential V app defined by 



(2.9) V ap Jx 



x when x G (— 7r\/2, 0), 
when x > 0. 



8 Airy/2' 
1 
2 

After the change of variables z = y/2x/(3ir) and the change of parameter k = 4/i/(37P\/3), 
we find the correspondence 

(2.10) - h 2 d 2 x + V app (x) ~ ^H toy (K)[z; d z ] + i 

where the toy model operator TAoy(ft) [z; d z ] is defined as: 



—z when z < 0, 



(2.11) U toy {K) = -K 2 d 2 z +W{z) with W(z 



This toy model invites us to recall the properties of the Airy operator. 



2.2.2. The Airy function and its zeros. Let us recall the basic properties of the Airy operator, 
i.e. the Dirichlet realization on L 2 (M_) of the operator — d\ — z. The electric potential tending 
to infinity when z — > — oo, this positive operator has compact resolvent. Thus, its spectrum 
can be described as an increasing sequence of eigenvalues tending to +oo. Let us use the 
traditional notation Ai for the Airy function. We recall that it satisfies: 

-Ai" + zA\ = 0. 

All along this paper, we will use A the reverse Airy function, i.e. A(z) = Ai(— z). We recall 
that A does not vanish on R_, is exponentially decreasing when z — > — oo and that its zeros 
(which are simple) form an increasing sequence of positive numbers tending to +oo. 

Notation 2.2. The n-th zero of A are denoted by z A (n). 
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If (A, is an eigenpair of the Airy operator, we have — ip'[ — zvpx = \ip\, hence the equation 
— ijj'l — (z + A)^a = 0. We deduce that there exists a number c(A) so that: 

^) = c(A)A(z + A). 

With those remarks, we can see that the spectrum of the Airy operator is {z/\{n),n > 1} and 
these eigenvalues are simple. 

2.2.3. Born-Oppenheimer approximation on the triangle. Finally, let us introduce the Dirich- 
let realization on L 2 ((— 7r\/2, 0)) of: 

(2.12) UboM = -h 2 d 2 x + - - . 

This operator is the Born-Oppenheimer "approximation" of the operator Cj r \{h) on the triangle 
Tri. The proof that it is actually an approximation will be done in Subsection 5.4 through 
the Feshbach projection: Indeed the operator 1-LQoj r \(h) has the same two-term eigenvalue 
asymptotics as the operator Cj r \(h) on the triangle. 

2.3. Asymptotic expansions of eigenvalues. We are now in position to state the results on 
eigenvalue expansion that we have proved in this paper. 

2.3.1. One-dimensional models. The lowest eigenvalues of the toy model (2.11) admit ana- 
lytic expansions with respect to k, 1 ^ 3 (when k is small enough): 

Theorem 2.3. For all N G N, there exists k > such that, for k G (0, kq), there exists 
at least N eigenvalues ofH toy (n) below 1. Denoting by \ y,n{^) the increasing sequence of 
these eigenvalues, we have the converging expansions for 1 < n < N and k small enough: 

+oo 

Atoy,n( K ) = ft 2 / 3 a.^ n K?fi with first coefficient a 0i n = z^n). 

3=0 

The corresponding eigenvectors have expansions in powers ofn 1 / 3 with the scales z/k 2 ^ 3 when 
z < and z/h when z > 0, see (3.7)-(3.8). 

As already mentioned we will meet in our investigation the Born-Oppenheimer approxima- 
tions of C Tr \(h) and £ Gui (/i). In order to compare the different asymptotics, let us state the 
result about the eigenvalues of V,BOjn(h): 

(2.13) A B o,Tri,n(/>) ~ V knh 2j/3 with /§ , n = I and {J 1>n = (4 7 rv / 2)~ 2/3 ^ A (n), 
and about the eigenvalues of W.BO,Gu\(h): 

(2.14) A B0 ,Gui,n(/i) ~ y^lj,nh j/3 with 7o , n = -, 7 ln = 0, and7 2 , n = (4irV2)- 2/3 z A (n). 

i>o 

Let us point out that this latter estimate will not be used to prove our main theorem (see 
Theorem 2.6) but somehow reflects that ?^BO,Gui(^) is an approximation of £c u i(^)- 
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2.3.2. Triangles. The lowest eigenvalues of the triangle Trig admit expansions at any order in 
powers of O 1 ^ 3 . We first state the result for the scaled operator Cj r \(h) introduced in §2.1.3: 

Theorem 2.4. The eigenvalues of Cj r \(h), denoted by \j r \,n{h), admit the expansions: 

\jn,n( h ) ~ Y]Pj,nh i/3 With 00,* = -, /?!,„ = 0, and fan = (^V2)- 2/3 Z A (n) , 
h— >0 ^— ' O 

i>o 

the terms of odd rank being zero for j < 8. The corresponding eigenvectors have expansions 
in powers ofh 1 ! 3 with the two scales x/h 2 / 3 and x/h, see (5.14). 

In terms of the physical domain Trig, we deduce immediately from the previous theorem 
that the eigenvalues of A^f , denoted by pj r -^ n (9), admit the expansions: 

A*Tri,»(0) ~ E^ 73 with = i 0i, n = 0, and /3 2 A n = 2(4vrv / 2)- 2/3 z A (n), 

j>0 

the coefficients /3^ n having the same properties as the /3j >n . Performing the dilatation: 

Xi = sin26'xi X2 = sin 2$ £2, 

we transform Trig into a new isosceles triangle with angle a = 29 and two sides with length c = 
2ix. Let us denote by p,j~. (a) its Dirichlet eigenvalues. It is easy to see that the eigenvalues 
satisfy the relation: 

0Tri,n(0) = (s in «)Vj>i,n(>)> 

so that we find back the result of [20, Theorem 1]. 

Remark 2.5. As it will be seen in the proof, the existence of a non-zero coefficient /3g n at the 
order 9 in the expansion of A-rri,n(^) reduces to the evaluation of an integral, see (5.12). If 
fan 7^ 0, there is a nonzero odd term after 0(a 2 ^ 3 ) in the asymptotics of /i^ 1 (a). 

2.3.3. Broken guides. Finally, we state the approximation result for the eigenvalues of the 
scaled operator £G u \(h) introduced in (2.5): 

Theorem 2.6. For all No, there exists ho > 0, such that for h G (0, ho) the No first eigenvalues 
of Cc u \{h) exist. These eigenvalues, denoted by \Q u \, n {h), admit the expansions: 

Acui.nO) ~ y^7j,n/i j/3 Withj ,n=Q, Jl,n = 0, and J 2 ,n = (47T\/2)~ 2/3 Z A (n) 
h— >0 ' ■* o 

i>o 

and the term of order h is not zero. The corresponding eigenvectors have expansions in powers 
ofh 1 ! 3 with the scale x / h when x > 0, and both scales x / h 2 ^ 3 and x / h when x < 0, see (6.12). 

Deducing the eigenvalues in the waveguide fig (Theorem 1.2) is an obvious consequence of 
this theorem. 
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2.4. Notation and terminology. The L? norm will always be denoted by || ■ ||, in general 
without mention of the integration domain. For a subset S C M. and a point pel, dist(S',p) 
is the distance between S and p, i.e. inf se5 \s — p\. 

We denote by &(A) the spectrum of a self-adjoint operator A, by <5 e ss(^4) its essential spec- 
trum, and by & d - ls (A) its discrete spectrum. An eigenmode (or eigenpair) of A is a pair (A, ip) 
with ip in the domain of A, such that Aip = Xip; then A is the eigenvalue and ip the eigenvector. 
A quasimode for A is a pair (A, tp) such that \\Aip — \ip\\ < £\\ip\\ with e small; A is the quasi- 
eigenvalue and ip the quasi-eigenvector. The spectral theorem implies that dist(@(A), A) < e. 



Domain 


Notation 


Variables 


Main operators 


Scaled Triangle 
Rectangle 
Half- strip 


Tri (2.7) 
Rec (5.2) 

Hst = M_ x (-1,1) 


(x,y) (2.2) 
(u,t) (5.1) 
(s,t) (5.5) 
(cr,t) (5.5) 


Crn(h) = -h 2 %-% 

C Rec (h) (5.3) 
Y.jtiJ' 2 ' " (5-7) 
Y.jMsjh* (5-8) 


Scaled half-guide 
Left half-strip 

Right half-strip 


n (2.3) 

HIef = M_ x (0,1) 
Hrig = R + x (0,1) 


(x,y) (2.2) 

(s,t) 

(a,t) 

(a,r) (6.2) 


cUh) = -h 2 dl-dl 

VjCrJr 1 ;! Notation 6.2 
Ei-A/"sf^ Notation 6.2 
Ei-^s?^ Notation 6.2 



Table 1 . Main notation for domains, variables and operators. 



3. Toy model in one dimension 

This subsection is devoted to the proof of Theorem 2.3 devoted to the spectral asymptotics 
of the operator "H t o y (^) defined in (2.11). This proof is divided into two steps. First, we 
construct quasimodes for %toy(ft), and second, we show that the lowest quasi-eigenvalues are 
the approximations of the lowest eigenvalues of "Htoy(^) of the same rank. 

3.1. Construction of quasimodes. In this section we prove in particular the following: 
Proposition 3.1. For all N e N*, there exists k > and C > such that for k E (0, k ): 

(3.1) dist(6 d i S (ft toy (/0), ^aH) < CVc, n = 1,---N . 

Proof. The basic tool for the proof is the construction of quasimodes and the application of the 
spectral theorem. Convenient quasimodes are given by power series in k 1 ^ 3 of profiles at the 
scales 

(3.2) s = k~ 2 ^z when z < (left) and a = k~ 1 z when z > (right). 



(3.3) Mz) 



(3.5) 
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More precisely we look for quasi-eigenfunctions ip K in the form: 

E, t/I'im, (•-)'••' :i when^<0 
Ej>o $ rig,i(^) K j/3 when z > , 
and quasi-eigenvalues in the form: 

(3.4) a K ~ k 2/3 a J Ki/3 as k — y 0. 

i>o 

The continuity conditions at z = provide the formal identities: 

«- a/8 E^o^wj(o)^ /s = ^E^o^ngjCo)^ 3 , 

and the formal eigen-equation is 

(3.6) -K 2 ^{z) + W{z)i> K (z)=a K i> K (z) zeR. 

• Determination of ceo. Collecting the terms in k 2//3 in (3.6) and using (3.3)-(3.5) we obtain: 
~K g ,o( a ) + $ ri g ,o(<r) = for a > 0, and $' rig0 (0) = 0, 

-*lef, ( S ) - S*|ef,o(s) = a *lef,o(a) for S < 0, and *|ef, O (0) = $rig, O (0). 

We deduce first that $ r ig,o = and thus ^i e f,o(0) = 0- This implies that a is a zero of the 
reverse Airy function A. At this stage we can choose a positive integer n, take a = z^{n) and 
\&ief,o as the corresponding normalized eigenfunction 

• Determination of ax. Collecting the terms in k, we get the equations: 

-^ g)1 + $ rig ,i = for a > 0, and $' rig)1 (0) = ^ ef>0 (0), 

-*;; f)1 - s*,ef,i - a *ief,i = ai^ief.o for s < 0, and *, e f,i(0) = $ rig ,i(0). 
We find first: 

$ rig ,i(<r) = -^ efi0 (0)e- CT . 

Moreover we obtain the existence of a number ol\ and of an exponentially decreasing \Pief.i 
solution of the second equation with the help of the following lemma: 

Lemma 3.2. Let n > 1. We denote by gc n ) an eigenvector of the operator —d 2 — s associated 
with the eigenvalue z&(n) and normalized in L 2 {R_). Let f = f(s) be a real function with 
exponential decay and let cGi Then there exists a unique aGS such that the problem: 

(-d 2 - s - z A (n)) g = f + ag {n) in R_, with g(0) = c, 

has a solution with exponential decay. There holds 

a = c ^(n)(°) - / f( s ) 9(n)(s) ds. 

• Further terms. A similar procedure can be reproduced at each step, providing the construc- 
tion of $ri g j» then and ^lefj, for any j > 2. 
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• Expressions for quasimodes. Relying on the previous iterative constructions we can set for 
all integer J > 

( J+2 

E*'ef,;(^)^ /3 when ,<0 



(3.7) ^\z 



j=0 
J+2 



^' /3 + ^eW +2 (0)^ /3 ^(9 when z > 



where x is a smooth cutoff function equal to 1 near 0. By construction, ^I J ' and its first 
derivative are continuous in z = 0. Moreover ipjf^ is exponentially decreasing as z — > ±oo. 
Therefore it belongs to the domain of T-L toy (k) . With this remark and taking the error introduced 
by x into account, we get for all k > 0: 

J+2 

(3.8) ^(n toy (K) - K 2/3 (z A (n) + J2 a J KJ/3 ))^l J] \\ <C{J,n,K )K l+J/ \ V/c < k . 

j'=i 

Hence 

||(^toy(K) - « 2/3 -2a(™)) ^k\\ < C(n,K ) K, Vk < Kq, 

and the spectral theorem applies. In particular, for k small enough, the discrete spectrum of 
H toy (K) is not empty since 6 e ss('Htoy(^)) = [1, +oo). □ 

Remark 3.3. We have proved in fact more than Proposition 3.1. The expression (3.7) of quasi- 
modes and corresponding estimates (3.8) will provide an asymptotic expansion for the eigen- 
vectors of Htoy(K), once one knows Proposition 3.4 below. 

3.2. Localization of the lowest eigenvalues. We now want to refine Proposition 3. 1 by prov- 
ing that the A toy ,n( K ) are power series with respect to k 1//3 and whose coefficients are given by 
(3.4). We begin to prove the following proposition: 

Proposition 3.4. For all No £ N*, there exists k,q > and C > such that for k £ (0, Kq): 

(3.9) lV,n(«) - k 2/3 *a Wl <Ck, n = 1, • • • N . 

Proof. Let N £ N*. As a consequence of Proposition 3.1, we have in particular that, for all 
k £ (0, k ), the first iVo eigenvalues A toy>n (fi;) (denoted by A n for shortness) exist and satisfy: 

(3.10) |A n | < C{N q )k 2 /\ «6(0,ko), n = l,-.. No- 
Let us denote by ip n an eigenfunction associated with A n so that (ip n , ip m ) = if n ^ m. For 
z<0we have: 

Thus, there exists a coefficient c n (n) such that: 

(3.11) Xpn(z) = C n (K)A(K~ 2/3 Z + K" 2/3 A n ), Z < 0. 

For z > we have the equation —K 2 tp^ = \ n ip n , hence the existence of d n (n) such that: 

(3.12) = d n (K)e~ K " zVI ^, z > 0. 
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The transmission conditions at z = imply: 

c n («;)A(«r 2 / 3 A n ) = d„(«), CnMK^A'fa- 2 / 3 ^) = -d n ( K )Vl- A n . 
This implies: 

(3.13) 
We infer: 



A(*T 2 / 3 A r 



A'(«T 2 / 3 A n ). 



|A(K" 2 / 3 A n )| < C7(iV ) k 1 / 3 . 

Since K~ 2 / 3 A n is bounded, see (3.10), and the zeros of the Airy function being isolated and 
simple, we deduce that for all n 6 {1, • • • , N Q }, there exists p = p(n, k) such that: 

|K- 2 / 3 A n -^ A (p)| < C7(iV )K 1 / 3 . 

Note that p is bounded too. It remains to prove that p = n for k small enough. In view 
of Proposition 3.1, it suffices now to prove than if k is small enough and n ^ m (with n, 
m < Nq), the integers p(n, k) and p(m, k) are distinct. Let us prove this by contradiction. 
Since the considered sets of integers n, m and p are finite, the negation of what we want to 
prove can be written as 

3m,n,p6N, V^i > 0, 3k 6 (0, k±) such that p(m, k) = p(n, k) — p. 

The eigenfunctions can be taken in the form: 

A(k~ 2/3 z + K" 2/3 A i ) when z < 



and we have 
A rescaling leads to 



for J = TYl Tl 

A(k- 2 / 3 A,) e -^'y/^ when z > , 
A(«T 2/3 z + K^ 2/3 A n )A(K" 2/3 z + «T 2/3 A m ) dz + 0(k 5/3 ) = 0. 



z<0 



[ A(z + *r 2/3 A n ) A(z + /sT 2/3 A m ) 
</2<o 



< C(N )k. 



By assumption, « 2 ^ 3 A n = z A (p) + (^(k 1 / 3 ) and k 2 ^ 3 X m = Za(p) + O^k 1 /' 3 ). For j = n, m, 
A being Lipschitz on (— oo, M] for all M, there exists D(N ) > such that for all z < 0: 

|A(z + k- 2/3 Aj-) - A(z + z A (p))| < J D(A r o)« 1/3 , for j = m,n, 

so that: 



'2<0 

We deduce: 



/ A(z + K- 2 / 3 X n )A(z + K-^ 3 X m )dz- [ A 2 (z + z A (p)) 

Jz<0 Jz<0 



dz 



< D(N )k 



1/3 



Vki > 0, 3k E (0, Ki) such that 



A 2 (z + z A (p))dz 



z<0 



< D(N )k 



A/3 



which leads to a contradiction and ends the proof of Proposition 3.4. 



□ 
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3.3. Proof of Theorem 2.3. Let us observe that Proposition 3.4 allows to separate the first N 
eigenvalues when k < k . Let us write 5 = k 1 ^. We let: 

so that A n (<5) is uniformly bounded for n = 1, . . . , No and 5 < k q . 

1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 

( 

Figure 4. The first two eigenvalues A toy ,i and A toy ,2 as functions of 5 = k 1//3 . 





0.48 0.49 0.5 0.51 0.52 0.53 0.54 0.55 0.56 0.57 0.58 



0.98 
0. 



Essential spectrum 
- First eigenvalue 
Analytic extension 



0.95 0.96 0.97 0.98 0.99 1 1.01 1.02 1.03 1.04 



Figure 5. The eigenvalues A toy ,i (left) and A toy ,2 (right) as functions of 5 
K 1 / 3 , zoom near the bottom of the essential spectrum. 



We rewrite (3.13) in the form: 
(3.14) A(A n (£)) 



S 



1 - S 2 K(S) 



A'(X n (5)). 
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We know that A is analytic and, using again the simplicity of its zeros, we can apply the 
analytic implicit function theorem near 5 = and for all n e {1, • • • ,N }, which, together 
with (3.7)-(3.8) and Proposition 3.4, ends the proof of Theorem 2.3. 

Remark 3.5. From (3.14), we can deduce that the X n (5) are solutions of the analytic equation: 
(3.15) (1-5 2 A)A(A) 2 -5 2 A'(A) 2 = 

This equation provides an analytic extension of the functions 5 A n (5), hence of A toyri = 
<5 2 A n (<5), in the sense of analytic curves. We represent in Figures 4 and 5 the first two eigenval- 
ues and their analytic extensions. Taking the continuity and monotonicity of the eigenvalues 
with respect to 5 into account, we can see that any branch which starts by 5 \(5) = 
5 2 z& + 0(S 3 ) represents an eigenvalue while \(S) is less that 1. Beyond 1, the Rayleigh quo- 
tient stays = 1, but the curve A (5) has an analytic extension as a continuation of a branch of 
roots of the equation (3.15). 

4. BORN-OPPENHEIMER APPROXIMATION FOR THE TRIANGLE 

This section is devoted to the analysis of V-BOjnih) defined in (2.12). We are going to 
prove: 

Theorem 4.1. The eigenvalues ofHsojr\{h), denoted by ABO,Tri,n(M> admit the expansions: 

ABO,Tri,nO) ~ J^Pj,nh 2j/3 , With Q ,n = I and $ 1>n = {AllV^^^Z^n). 

i>o 

Again, the proof is essentially organized in two steps. The first step is the construction of 
quasimodes which proves that quasi-eigenvalues are close to true eigenvalues. The second step 
uses Agmon type exponential localization for true eigenvectors to prove that true eigenvalues 
are close to quasi-eigenvalues. 

4. 1 . Quasimodes. In this subsection, we construct quasimodes and prove the proposition: 
Proposition 4.2. For all N e N*, there exists h > and C > such that for h G (0, ho): 

(4.1) dist(6 dis (HBo,TriW)^ + ^ /3 (4vrv / 2)- 2/3 z A (^)) < Ctfl\ n= 1,---N . 

Proof. The proper scale in x is h 2 ^ 3 as can be seen by approximating the potential in x = 
by its tangent and recognizing the Airy operator. Thus, we will construct quasimodes ^ as 
functions of s = h~ 2 l 3 x: We look for quasimodes (A^, iph) m the form of series 

in order to solve HBOjriifyiph — Xh*l>h in the sense of formal series. A Taylor expansion at 
x = of the potential V yields: 

HnoMh) - ~ h2d l + J2 with v ° = \ and y i = -ir-m> 
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which, in s variable, becomes 

(4.2) n B0J n(h) ~ i + /> 2 / 3 ( - <9 S 2 + ^a) + ]T // 2 -' 

i>2 

The construction of the terms (3j and ^ is similar (even simpler) than for Proposition 3.1. 

• The expansion (4.2) yields that /3 = |, and collecting the terms in h 2 /' 3 and we obtain: 

(4.3) _^'( 5 )__i_^ ( s ) = ^ ( s ) Vs<0 and tf (0) = 0. 

Thus for any chosen positive integer n we can take (3i = (47tv / 2) _2 ^ 3 ^a(^) together with 

(4.4) * (a) = A((4vr V / 2)~ 1/3 s + z A {n)). 

• Collecting the terms in /i 4 / 3 we obtain 

-^ , 1 / (s) + V / is^ 1 (s)-/3i^i(s) = /3 2 ^o-V A 2 s 2 ^o Vs<0 and ^(0) = 0. 

The compatibility condition states that ^(^Oj ^o) = ^(s^cb ^o)- This determines (3 2 and 
implies the existence of a unique solution £ L 2 (IR_) such that \l/ ) = 0. 

• This procedure can be continued at any order and determines 0j, *ffj) at each step. This 
construction depends on the choice of the integer n and can be done for any positive integer n. 

• To conclude, we consider a smooth non-negative cutoff function x' ef satisfying: 

(4.5) x ]ef (x) = 1 for x £ ( - -j=, +oo) and ^ e \x) = for x < -it, 

and introduce for any J > the quasimode (/3[ J ', ^1^) with: 

(4.6) ,/ ir -hir^ and ^(x) = x' ef (x) £ *i (t^)^ 3 - 

3=0 3=0 

Thanks to this cut-off -0^ satisfies Dirichlet condition in —tx^/2, and in by construction. 
Using the exponential decay of x t-> ^ j{h~ 2 ^x) and the definition of \T/j and we get for 
any h > the existence of C(n, J, h ) > such that: 



(4.7) 



(-HbojM - /f V1 J1 < C(n, J, h ) /i 2(J+1)/3 , V/i £ (0, ho) 



This proves the existence of quasimodes at any order and ends the proof of Proposition 4.2. □ 

4.2. Agmon estimates. In this subsection, we prove Agmon estimates (see [1, 2]) for the 
eigenfunctions of Heoj^h). The role of Agmon estimates is to replace an explicit knowledge 
of the solution at infinity like in (3.1 1)-(3.12) by suboptimal exponential estimates. 

Here we prove two kinds of estimates: near x = and near x = —n\^2. In the analysis of 
triangles (cf. Section 5.2), we will meet the same estimates. Let us consider an eigenpair (A, ip) 
of %BO,Tri(^)- The Agmon identity writes, for some Lipschitz function $ to be determined: 

(4.8) f h 2 \d x {e*^)\ 2 + V{x)\e*tP\ 2 -h 2 \<S>'e*\ 2 -\\{e*4j)\ 2 dx = 0. 
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It is a consequence of Proposition 4.2 that the lowest N eigenvalues A of "HBO,Tri(^) satisfy: 
(4.9) |A-||<r /i 2 / 3 , 
for some positive constant T depending on N . 



4.2. 1 . Agmon estimates near x — 0. We use (4.8) and the convexity of V to get the inequality: 

-o 

-ttV2 ' XK V 8 47T\/2 / 



h 2 \d x (e*if>)\ 2 + ( i - — ^= ) |e^| 2 - /i 2 |$'e*| 2 - \\{e*i/>)\ 2 dx < 0. 



f %|e*Vf - fc 2 |$'e*| 2 - C/z 2 / 3 |(e*V)| 2 ^ < 0. 

J-nV2 47TV2 



With (4.9), we deduce: 

»o 

This leads to the choice 

for a number 77 > to be chosen small enough. We get: 







\x\ 9 2 



-7] \X\ 



Ch 2 ^ \e^\ 2 dx < 0. 



-*VS \47r v / 2 4 
For 77 small enough, we obtain the existence of fj > such that: 



/ (77 |a; I -C7i 2/3 ) |e*^| 2 dx < 0. 



Splitting the integral into the parts < x < —Dh 2 ^ 3 (where $ is unbounded) and 

—Dh 2 ^ 3 < x < (where $ is bounded) with f\D — C = d > 0, we find: 



-Dh 2 / 3 p-Dh 2 ' 3 



/—un. i" r—un i 

d0 3 \e*ip\ 2 dx< / (fj\x\ -Ch 2 / 3 ) \e*i/)\ 2 dx 



tV2 J 

< I (fj\x\ + Ch 2 ' 3 ) |e*^| 2 dx < Ch 2 ' 3 I |z/f dx. 

J-Dh 2 / 3 J-Dh 2 / 3 

We deduce the proposition: 

Proposition 4.3. Let T > 0. There exist h > 0, C > and r/ > swc/z that for h G (0, /i ) 
and all eigenpair (A, ^) ofH&oj r \{h) satisfying |A — || < T h 2 / 3 , we have: 

e mh~^ ( |^|2 + \0 3 d x ^\ 2 ) dx < Co ||^||2_ 
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4.2.2. Agmon estimates near x = —n y/2. We use again (4.8) and (4.9): 

f h 2 \d x {e^)\ 2 + ( ^ - \) |eV| a - h 2 \&e^\ 2 - C0 3 \(e*i/>)\ 2 dx < 0. 
J-nV2 V4(x + vrv / 2) 2 8/ 

We take: 

$(z) = -phT 1 In (D~ 1 (x + nV2)), 
where we choose p E (0, |) so that there holds: 

f 2 ((t - p2 ) (x + 7rv ^~ a " |e ^' 2 - ^ a/3 l( e< W dx ^ °. 

and D > large enough so that 

(T-')*-i> tt 

Then we split the integral into the parts —tt^/2 < x < — 7r\/2 + D (where $ is unbounded) 
and —ix\[2 + D < x < (where $ is bounded) and the same procedure as in the previous 
paragraph leads to the proposition: 

Proposition 4.4. Let T > and po £ (0, §). TTiere exxst /io > 0, Co > smc/i that for any 
h E (0, h ) and all eigenpair (A, ip) ofT-Lsojnih) satisfying |A — ~| < T h 2 ^ 3 , we have: 

o , x 

(x + 7tv / 2)" Po//i ( \ip\ 2 + I^Vf) dx < CoU\\ 2 . 

4.3. Proof of Theorem 4.1. Let us fix No and consider the N first eigenvalues of "Hbo,tn(^) 
denoted by \ n = \BOjn,n(h). For each n E {1, • • • iV }, we choose a normalized ipn in the 
eigenspace of A n so that (ip n , ip m ) = for n ^ m. Let us introduce the space: 

&N (h) = span (-01, • • • , ipN )- 
We recall that, for h small enough, (4.9) holds. We can write: 

%BO,Tri(fr)0n = A n ^ n 

so that (the ip n are orthogonal in L 2 and for the quadratic form), for all ip E (£jv (/i): 

QBO,Tr\,h(i>) < ^Afoll^ll 2 - 

For e small enough we introduce a smooth cutoff function \ being for \x + ix\pl\ < e and 
1 for \x + tx\[2\ > 2eq. Proposition 4.4 implies that: 

Qbo ,Tri,h 

Then, the convexity of the potential allows to write: 

-h 2 d 2 x - ^= x + < (Atvo + 0(/O)«ir, 

where we have used the convexity. The dimension of x&N {h) ls s ° that, with the properties 
of the Airy operator and the mini-max principle, we get: 

I + (4vrv / 2)- 2/3 2A (iVo) < X No + 0{h°°). 



20 



MONIQUE DAUGE AND NICOLAS RAYMOND 



This is true for all fixed N and, combined with Proposition 4.2, provides the separation of the 
lowest eigenvalues of V,Bo,Jr\(h): 



Finally, with (4.6)-(4.7), we obtain Theorem 4.1. 

4.4. Born-Oppenheimer approximation for the waveguide. Let us end this section by an 
informal analysis of the spectrum of the operator "HBO,Gui(^) defined in (2.8). This investiga- 
tion is not necessary in our way to prove Theorem 2.6, but it already gives a flavor of the ideas 
to analyze CG u \(h). We can obtain the asymptotic expansions of the eigenvalues of V,BO,Gu\(h) 
by combining the analysis for % t0 y{h) and for %bo,Tn(^)- Indeed we can perform a quasimode 
construction like for H t oy(h) and 7^BO,Tri(^) by solving a transmission problem between the 
negative half-axis and the positive half-axis. For that purpose, we establish the following Ag- 
mon type estimate which states that the eigenfunctions of W.BO,Gu\(h) do not penetrate in the 
region x > more than at the scale h. 

Proposition 4.5. Let (A, -0) be an eigenpair ofHBO,Gu\(h) such that |A — || < Ch 2 ^ 3 . There 
exist a > 0, ho > and C > such that for all h G (0, ho), we have: 



Thanks to the latter estimate we cut off the part of the eigenfunctions living on x > modulo 
a remainder of order 0(h°°). This allows the comparison of %BO,Gui(^) with Hbo,th(^) an d 
provides the proof that all the lowest eigenvalues of ^BO,Gui(^) are described by the quasimode 
construction. In Section 6.2 a similar analysis will be done for the whole waveguide. 



The aim of this section is to prove Theorem 2.4. As usual, the proof will be divided into two 
main steps: a construction of quasimodes and the use of the true eigenfunctions of Cj r ;(h) as 
quasimodes for the Born-Oppenheimer approximation in order to obtain a lower bound for the 
true eigenvalues. 

We first perform a change of variables to transform the triangle into a rectangle: 





5. Triangle with Dirichlet boundary condition 



(5.1) 



u = x G (-7iV2, 0), t 




y 



1,1). 



so that Tri is transformed into 



(5.2) 



R ec = (-7iV2,0) x (-1,1). 



The operator Cj r \(h) becomes: 



(5.3) £ Rec (h)(u,t;d u ,d t ) = -h 2 (d u - 



t 




1 




U + 7T\/2 
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with Dirichlet boundary conditions on dRec. The equation Cy r \(h)t/j h = finish is transformed 
into the equation 

jC Rec (h)ip h = (3 h ip h with iph(u,t) = ip h (x,y). 

5.1. Quasimodes. This subsection is devoted to the proof of the following proposition. 

Proposition 5.1. There are sequences (/3j,n)j>o far any integer n > 1 so that there holds: 
For all JVo Gl and J G N, there exists ho > and C > such that for h G (0, ho) 

j 

(5.4) dist(6 dis (£ Tri (/*)), Ys^ h3 ' Z ) ^ ChiJ+m i n = l,---N . 

j=0 

Moreover, we have: 0o >n = |, 0i j7l = 0, and /3 2 ,n = (^\/2)~ 2 / 3 z&(n). 

Proof. We want to construct quasimodes ((5 h , ip h ) for the operator Cj r \(h)(d x , d y ). It will be 
more convenient to work on the rectangle Rec with the operator C Rec (h)(u,t;d u ,d t ). We 
introduce the new scales 

(5.5) s = h~ 2/3 u and a = h~ l u, 
and we look quasimodes ((3h,iph) in the form of series 

(5.6) (3 h ^Y,^ ylZ and MM)~^(^M) + ^M))^ /3 

in order to solve C Rec (h)tjjh = Ph^h in the sense of formal series. As will be seen hereafter, an 
Ansatz containing the scale h~ 2 l 3 u alone (like for the Born-Oppenheimer operator 'Hbo,™^)) 
is not sufficient to construct quasimodes for C Rec (h) . Expanding the operator in powers of h 2 / 3 , 
we obtain the formal series: 

(5.7) £ Rec (h)(h 2/3 s, t; h~ 2/3 d s , d t ) ~ V C 2j h 2j/3 with leading term C = --^d t 2 

j>o 

and in powers of h: 

(5.8) £ Rec (h) (ha, t; h~ l d a , d t ) ~ V Af 3j h j with leading term Af = -d 2 a - Ar<9 2 . 

In what follows, in order to finally ensure the Dirichlet conditions on the triangle Tri, we will 
require for our Ansatz the boundary conditions, for any j G N: 

(5.9) * J -(0,t) + $ i (0,t) = 0, — 1 < t < 1 

(5.10) * i (s,±l)=0, s<0 and $ i (a,±l) = 0, a < 0. 

More specifically, we are interested in the ground energy A = | of the Dirichlet problem for 
£ on the interval (—1,1). Thus we have to solve Dirichlet problems for the operators Mq — \ 
and jC — | on the half-strip 

(5.11) Hst = R_ x (-1,1), 
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and look for exponentially decreasing solutions. The situation is similar to that encountered in 
thin structure asymptotics with Neumann boundary conditions. The following lemma shares 
common features with the Saint- Venant principle, see for example [13, §2]. 

Lemma 5.2. We denote the first normalized eigenvector of X on ((— 1, 1)) by c : 

c (t) = cos f^j . 

Let F = F(a,t) be a function in L 2 (Hst) with exponential decay with respect to a and let 
G G iJ 3 / 2 ((— 1, 1)) be a function oft with C7(±l) = 0. Then there exists a unique 7GI such 
that the problem 

Af - - ) $ = F in Hst, $0,±1) = 0, $(0,t) = G(t) +7C (t), 



admits a (unique) solution in H 2 (Hst) with exponential decay. There holds 

7 



F(a,t)aco(t)dadt- [ G{t)c {t)dt. 
1 J -1 



The following two lemmas are consequences of the Fredholm alternative. 

Lemma 5.3. Let F = F(s, t) be a function in L 2 (Hst) with exponential decay with respect to 
s. Then, there exist solutionis) \I/ such that: 

Co-- )* = F in Hst, tf(s,±l) = 



if and only if (F(s, •), cq) = for all s < 0. In this case, ^f(s,t) = ^ ± (s,t) + g(s)co(t) 
where ^Z" 1 satisfies (^(s,-),^). =0 and has also an exponential decay. 



Then, we will also need a rescaled version of Lemma 3.2. 

Lemma 5.4. Let n > 1. We recall that zj\{n) is the n-th zero of the reverse Airy function, and 
we denote by 

g (n) = A((47rv / 2)" 1/3 s + z A (n)) 

the eigenvector of the operator — <9 2 — (47rv / 2)~ 1 s with Dirichlet condition on EL_ associated 
with the eigenvalue (47r\/2)~ 2 ^ 3 z&(n). Let f = f(s) be a function in L 2 (1L) with exponential 
decay and let c6i Then there exists a unique [3 G M such that the problem: 



(47rv / 2)- 2/3 2 A (^)^) 9 = f + (3g { n) in R_, with g{0) = c, 



has a solution in if 2 (R_) with exponential decay. 

Now we can start the construction of the terms of our Ansatz (5.6). 
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• Terms in h°. The equations provided by the constant terms are: 

£ ^o = A)^o(s, t), A/> = A)$o(s, *) 

with boundary conditions (5.9)-(5.10) for j = 0, so that we choose f3 = | and W (s,t) = 
g (s)c (t). The boundary condition (5.9) provides: $ (0, £) = — go(0)c (t) so that, with 
Lemma 5.2, we get go(0) = and $ = 0. The function go(s) will be determined later. 

• Terms in h 1 ^ 3 . Collecting the terms of order h 1 / 3 , we are led to: 

(£ - A))*i = ft*o - £1*1 = /9i*o, (A/"o - A))$i = /9i$o - M$i = 

with boundary conditions (5. 9)-(5. 10) for j = 1. Using Lemma 5.3, we find /?i = 0, \l/i(s, t) = 
^(s)co(t),^(0) =0and$! = 0. 

• Terms in h 2 ^ 3 . We get: 

(£ - /3 )^2 = /? 2 *o - £2^0, (A/o - A,)$ 2 = 0, 

where £ 2 = — ^ + <9 t 2 and with boundary conditions (5.9)-(5. 10) for j = 2. Lemma 5.3 
provides the equation in s variable 

((f3 2 ^o-C 2 ^o(s,-)),c ) t = 0, s<0. 

Taking the formula ^J/q = 9o(s)co(t) into account this becomes 



This equation leads to take /3 2 = (47r\/2)~ 2 / 3 2:A(n) and for g the corresponding eigenfunction 
gf( n ). We deduce (£ - A))*2 = 0, then get * 2 (s, t) = #2(»co(£) with g 2 (0) = and $ 2 = 0. 

• Terms in h 3 ^ 3 . We get: 

{C - A))^ 3 = A^o + /9 2 *i - £2*1, (A^o - /?o)$3 = 0, 

with boundary conditions (5.9)-(5.10) for j = 3. The scalar product with c (Lemma 5.3) 
and then the scalar product with g (Lemma 5.4) provide f3 3 = and g x — 0. We deduce: 

*s(s, *) = ^s(s)cb(*), and #,(0) = 0, $ 3 = 0. 

• Terms in /1 4 / 3 . We get: 

(£ - /3 )^4 = /3 4 *o + /?2*2 - £4*0 - £2*2, (A/" - /3 )$4 = 0, 

where 



a 3 2o2 



and with boundary conditions (5.9)-(5. 10) for j = 4. The scalar product with Co provides 
an equation for g 2 and the scalar product with g determines /3 4 . By Lemma 5.3 this step 
determines ^4 = \&f + c (t)g4(s) with a non-zero ^4 and #4(0) = 0. Since by construction 
(^i(0, •), c ). = 0, Lemma 5.2 yields a solution $ 4 with exponential decay. Note that it also 
satisfies ($4(0-, •), c ) f = for all a < 0. 
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• Terms in h 5 ^ 3 . We get: 

{C - /3 )*5 = + /3 2 ^s - C 2 V 3 , (Mo - /3 )$5 = 0, 

and with boundary conditions (5.9)-(5.10) for j = 5. We find /3 5 =0,g 3 = 0, ty 5 = g 5 (s)co(t), 
g 5 {0) = 0, $ 5 = 0. 

• Terms in h 6 ^ 3 . We get: 

(Co - A))^6 = A>^0 + /?4^2 + /?2^4 - C 2 ^4 - £4^2, CM) - /3 )$6 = /?2$4, 

and with boundary conditions (5.9)-(5.10) for j = 6. This determines (3 e , g$, \l/ 6 = \I/g + 
co(*)^6(s), #6(0) = 0, and $6 with exponential decay due to the orthogonality of $4 to c . 

• Terms in h 7 ^ 3 . We get: 

(£ - A>)*7 = £7*0 + £2*5 - £2*5, (A/" - /3 )$7 = -A/" 3 $4, 

where 



2 .„ „ a 

o 

t 5 



and with boundary conditions (5. 9)-(5. 10) for j = 7. We take f3 7 = 0, g 5 = 0, ^7 = g 7 (s)c (t). 
Then, Lemma 5.2 induces a value for the trace #7(0) so that there exists $ 7 with an exponential 
decay. Note that if there holds: 

(5.12) / (M^ A )(a,t)aco{t)dadt ^ 0, 

JHst 

we would deduce by Lemma 5.2 that #7(0) 7^ 0. 

• Terms in h 8 ^ 3 . We get: 

(£ - A,)*8 = fo*0 + /3 6 *2 + /?4^4 + #>*6 - £ 8 *0 ~ £ 6 *2 ~ £ 4 *4 ~ £ 2 ^6, 
CM) - A))$8 = /?4$4 + /3 2 $6- 

This determines /3 8 , g e and ty$ = + c ofi , 8, the trace gs(0) and the exponentially decreasing 
solution $ 8 . 

• Terms in h 9 ^ 3 . We get: 

(£ - /3 )^9 = A^O + /?2*7 - £2^7, (M) - A))$9 = /?2$7 ~ .M»$6- 

We find /3 9 , # 7 and then * 9 = ^ + c g 9 and 5(9(0), $ 9 . Note that if g 7 (0) ^ 0, i.e. if (5.12) 
holds, we would deduce that 0g 7^ 0. 

• Continuation. The construction of the further terms goes on along the same lines. This 
leads to define the quasimodes for Cj r \(h): 

j 



(5.13) Pr = Y,Pj } * 

3=0 

( , 14) ^^ {x) ±f H ^^ hHl ^ 
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where % lef is defined in (4.5). The conclusion follows from the spectral theorem. □ 

5.2. Agmon estimates. On our way to prove Theorem 2.4, we now state Agmon estimates 
like for V.BOjr\(h). Let us first notice that, due to Proposition 5.1, the lowest eigenvalues of 
Cj r \(h) still satisfy an estimate like (4.9). It turns out that we have the following lower bound, 
for all ip G Dom(Q T ri,/i): 

QjnAi')> [ h 2 \d x i;\ 2 + - U\ 2 dxdy. 

An 4(z + 7TV2) 2 

Thus, the analysis giving Propositions 4.3 and 4.4 applies exactly in the same way and we 
obtain: 

Proposition 5.5. Let T Q > 0. There exist h > 0, Co > and r] Q > such that for h G (0, ho) 
and all eigenpair (A, iji) ofCj ri (h) satisfying |A — || < T h 2 ^ 3 , we have: 

J e voh-i\x\*/* + |/,2/3^|2j dxdy < C ||^||2_ 

Proposition 5.6. Let T > 0. There exist h > 0, C > and p > swc/z that for h G (0, /i ) 
and a// eigenpair (A, ^) ofCj r \(h) satisfying |A — || < r /i 2 / 3 , we &ave; 

y (x + 7rV2)-^(|V| 2 + |/i^| 2 ) ^dy < C \m 2 . 

5.3. Approximation of the first eigenfunctions by tensor products. In this subsection, we 
will work with the operator £r 6C (/z) rather than £j r \ (h) . Let us consider the first N eigenvalues 
of £.R ec (h) (shortly denoted by A n ). In each corresponding eigenspace, we choose a normalized 
eigenfunction ip n so that (?/v, $ m ) = if n ^ m. As in Section 4.3, we introduce: 

&N (h) = span(^i, . . .,i>N )- 
Let us define Q^ ec the following quadratic form: 

QRecW = J (^2 l^l 2 " ^l 2 ) ( U + ^ dudt, 

associated with the operator £p ec = Id u g) (—^d 2 — |) on L 2 (Rec, (u + n\/2)dudt). We 
consider the projection on the eigenspace associated with the eigenvalue of —j^d 2 — |: 

(5.15) n $(u,t) = {4>(u,-),co) t co(t), 

where we recall that c (t) = cos . We can now state a first approximation result: 

Proposition 5.7. There exist h > and C > such that for h G (0, /i ) a// ip E <£n (h): 

o < QlM < Ch 2 ^\\ 2 

and 

||(id - n )^|| + ||<9 t (id - n )^|| < (7^/3^11. 

Moreover, Yl : <£N (h) —> Ho(<E]y (h)) is an isomorphism. 
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Proof. If ijj = tjj n , we have: 

QRec,h(Ai) = A n ||^„|| 2 . 

From this we infer: 

l|2 



The orthogonality of the 'ip n (in L 2 and for the quadratic form) allows to extend this inequality 

to^ £ <S No (h): 

This clearly implies: 

QLW < ct* 2/3 IHI 2 - 

n ^ being in the kernel of £r 6C , we have: 

Q RecW = Q Rec((W-n )^). 

If we denote by ^2 the second eigenvalue of the ID operator —-^d 2 ~ \-> we § et by the 
min-max principle: 

Q Rec ((id-n )^)>/i 2 ||(id-n )^|| 2 . 

Now the conclusions are standard. □ 

5.4. Reduction to the Born-Oppenheimer approximation. In this section, we prove The- 
orem 2.4 by using the projections of the true eigenfunctions (H ip n ) as test functions for the 
Born-Oppenheimer approximation. Let us consider an eigenpair (A, tp) of Cj r \{h) such that 
(4.9) holds. We let ip(u, t) = ip(x, y). Then, (A, ip) satisfies: 

2 2td u d t 2td t t 2 d 2 \ 1 1 ; . ; 

, = H = I =— w =— <x w = Aw. 

u + {u + Txy/2) 2 {u + ivV2) 2 J (u + ivV2) 2 

The main idea is to determine the (differential) equation satisfied by n ?/>. In other words we 
will compute and control the commutator between the operator and the projection n . For 
that purpose, a few lemmas will be necessary. The first one is an estimate established in the 
original coordinates (x, y) in the triangle Tri: 

Lemma 5.8. For all k £ N, there exist h > and C > such that we have, for h £ (0, ho): 

[ (x + irV2y k \d y iP\ 2 dxdy < C\\^\\ 2 . 

JTr\ 

Proof. The equation satisfied by ip is: 

{-h 2 dl - d 2 )^ = 

Multiplying by (x + n\/2)~ k , taking the scalar product with ip and integrating by parts we find: 

/ (x + nV2)- k \d y iP\ 2 dxdy < C [ (x + 7r V / 2)" fe (|Vf + h 2 (x + ti^2)~ 1 \^\ \d x ^\) dxdy. 
Using the Agmon estimates of Proposition 5.6 with po/h> fc + lwe deduce the lemma. □ 
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We can now prove: 

Lemma 5.9. There exist h > and C > such that we have, for h G (0, ho) 
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[u + tx\/2) Hdudttjj, co(t) 



< Ch 



-i 



L 2 {du) 



Proof. Integrating by parts in t for any fixed u G (— 0), we find: 

{{u + TxV2)'Hd u d t ^, c (t))J <C j {u + 7iV2)~ 1 \d J\dt 



1 \ 1/2 

2|o ?i2 



<C / (« + 7rV2)- 2 |9 u Vr^ 



To have the lemma, it remains to prove that 



/ (u + ixV2Y 2 \d u i)\ 2 dudt < Ch~ 2 [ 

JRec JRec 



We have: 



/ (u + nV2)- 2 \d u 4j\ 2 dudt = / {x + 7iV2)- 3 (d x + 



X 



dudt. 



dxdy 



and we apply Lemma 5.8 to control the term in d y . We end the proof using the Agmon esti- 
mates of Proposition 5.6. □ 

The same kind of computations yields: 
Lemma 5.10. There exist h > and C > such that we have, for h e (0, ho): 



(u + 7rV2)- 2 td t ip, c (t) 



L 2 (du) 



<cu\\- 



Finally, we have: 

Lemma 5.11. There exist ho > and C > such that we have, for h G (0, ho): 



(u + nV2y 2 t 2 dt^, c (t) 



L 2 (du) 



< c 



From Lemmas 5.9, 5.10 and 5.11, and from Proposition 5.7, we infer: 

Proposition 5.12. Let T > 0. There exist h > and C > such that for h G (0, h ) and all 

eigenpair (A, ip) of Cj r \{h) satisfying |A — \\ < T h 2 / 3 , we have: 



2q2 



-h 2 d, 



U 4(li + TTV^) 2 



- A U ip 



<Ch\\U ^\\. 
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• Proof of Theorem 2.4. We deduce, from Proposition 5.12, for all n 6 {1, 



No}: 



h 2 dl + 



7T 



U tjj r 



" 4(u + vr^) 2 
From this inequality, we infer, for all tp E <£jv (/i): 



< (XTn, NQ (h) + Ch)\\U ^ n \\. 



-h 2 d 2 



7l~ 



n ^ 



< (\jn,N (h) + ch)\\n ^\\ 



u 4(u + 7rV2) 2 . 
and thus: 

Qbo ,Tn,h (n V) < (Xrn, No (h) + Ch)\\U ^\\. 
It remains to apply the min-max principle to the iVo dimensional space U. (£.N (h) (see Propo- 
sition 5.7) and Theorem 4.1 to get the separation of eigenvalues. Then, the conclusion follows 
from Proposition 5.1. 



6. ElGENPAIR ASYMPTOTICS FOR THE WAVEGUIDE 

In this section, we prove Theorem 2.6. Firstly, we construct quasimodes and secondly we 
use Agmon estimates reduce to the triangle case. On the left, £ Gu \(h) writes, in the coordinates 
(u, t) defined in (5.1): 



(6.1) 



C'g ] {h) = -h 2 [d u - 



II 



a - 



1 r B] 
(u + ttV2) 2 



U — X, T 



y-x 



7l 



V2 



i 



and on the right, we let: 
(6.2) 

and the operator writes: 

The integration domain is {—iwfl, +oo) x (0, 1) = Cl\ ef U f2 rig where: 

Qief = (-ttV2, 0) x (0, 1) and fi rig = (0, +oo) x (0, 1). 

The boundary conditions are Dirichlet on (0, oo) x {0} U (— 7r\/2, oo) x {1} and Neumann on 
(-7iV2,0) x {0}. 

6.1. Quasimodes. The aim of this subsection is to prove the following proposition: 

Proposition 6.1. For any n > 1, there exists a sequence (7j, n ) such that, for all No e N and 
J G N, there exists h > and C > such that for h e (0, ho): 



j 

(6.4) dist(© dis (£ Gui (/>)), J2^ hj/3 ) ^ Ch(J+1)/ 

3=0 

Moreover, we have: 7o, n = §, Ji,n = and^2,n = (47tv / 2) _2 ^ 3 -2a 



n = l,---N n . 



n , 
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6.1.1. Preliminaries. 

• Ansatz, boundary and transmission conditions. In order to construct quasimodes for £q u \ {h) 
of the form (7^, ip h ), we use the coordinates (u, t) on the left and (u, r) on the right and look 
for quasimodes iph(u, t, r) = iph(x, y). Such quasimodes will have the form on the left: 

(6.5) Mu,t) ~ kj/3 (*lefj(^" 2/3 «,<) + ^efJ^M)) , 

i>o 

and on the right: 

(6.6) i/; rig (u, t)~J2 h i '^ s j{br 1 u, t) 
associated with quasi-eigenvalues: 

We will denote s = h~ 2 ^ 3 u and a = h^u. Since iph has no jump across the line x = 0, we 
find that r tp\ ef and r tp r]g should satisfy two transmission conditions on the line u = 0: 



Vte(0,t) = -0ri g (O,O and (<9„ ^<9, 



) ^lef(0,t) = (0 u --^=)^ rig (O,t), 



for all t E (0, 1). For the Ansatze (6.5)-(6.6) these conditions write for all j > 

(6.7) *lefj(0,t) + $|efj(0,t) = $rigj(0,t) 

(6.8) 9 CT $| efj (0,t) + ^^(O,*) j=$| rfJ _ 3 (0,t) ^W lef)i _ 3 (0,t) 

7TV^ 7TV 2 

= 9 CT $ rig , i (0,t)-^=$ rig , i _ 3 (0,t), 

7T V Z 

where we understand that the terms associated with a negative index are 0. 

Notation 6.2. We still set s = h~ 2 / 3 u and a = h~ 1 u. Like in the case of the triangle Tri, the 
operators C ] q u] and ^Gui' c f- (6.1)-(6.3), written in variables (s, t) and (a, t) expand in powers 
of h 2 / 3 and h, respectively. Now we have three operator series: 

• £jsfjj(/i)(/i 2 / 3 s, £; h~ 2 / 3 d s , d t ) ~ Z]j>o £-2jh 2j / 3 . The operators are the same as for Tri, 
but they are defined now on the half-strip Hlef := (—00, 0) x (0, 1). 

• £^{h){ha, t; h~ l d a , d t ) ~ E^o^f^' defined on Hlef - 

. cgMiha^h-^dr) -E,> K gfl3 defined on Hri S := (°'°°) x (Q' 1 )- 

We agree to incorporate the boundary conditions on the horizontal sides of Hlef in the defini- 
tion of the operators Cj, A/J ef , and M" g : 

• Neumann-Dirichlet d t ^(s, 0) = and *(s, 1) = (s < 0) for C jt 

• Neumann-Dirichlet d t <f>(a, 0) = and V(a, 1) = (a < 0) for A/"J ef , 

• Pure Dirichlet 0) = and (a, 1) = (a > 0) for A/J ig . 
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Note that 



< ef = - d "~^ d " and K^-^-^a 



• Dirichlet-to-Neumann operators. Here we introduce the Dirichlet-to-Neumann operators 
T ng and T lef which we use to solve the problems in the variables (a, t). We denote by / the 
interface {0} x (0, 1) between Hrig and Hlef. 

On the right, and with Notation 6.2, we consider the problem: 

(K S " ^ $ri g = in Hrig and $ rig (0,t) = G(t) 

where G G H^ 2 (I). Since the first eigenvalue of the transverse part of A/^ rig — | is positive, 
this problem has a unique exponentially decreasing solution $ rig . Its exterior normal derivative 
—da&ng on the line I is well defined in H~ l ' 2 {I). We define: 

T" g G = d $ • = -d $ ■ 

We have: 

(T^G,G)^Q r]g ^ rig )>C\\G\\" 
On the left, we consider the problem: 

(No f - -)$,ef = in Hlef and $, ef (0,t) = G(t) 
where G G H$(l). 

1 /2 

For all G G i/ g (/) such that n C7 = (where n is defined in (5.15)), this problem has a 
unique exponentially decreasing solution $| ef . Its exterior normal derivative d a &\ e f on the line 
/ is well defined in H~ l l 2 {I). We define: 

We have: 

(T lef G,G> = Q lef ($ lef )>0. 

Proposition 6.3. The operator T ng + T lef LTi is coercive on H^ 2 (I) with ni = Id - LT . In 

particular, it is invertible from Hqq 2 (I) onto H~ l l 2 (I). 

This proposition allows to prove the following lemma which is in the same spirit as Lemma 
5.2, but now for transmission problems on Hlef U Hrig (we recall that c (t) = cos(|t)): 

Lemma 6.4. Let F\ ef = F\ e f(a,t) and F rig = F rig (o", r) be real functions defined on Hlef 

and Hrig, respectively, with exponential decay with respect to a. Let G° G Hqq 2 (I) and 
H G H^ 1 / 2 ^) be data on the interface I = <9Hlef H <9Hrig. Then there exists a unique 

1 /2 

coefficient C, G M. and a unique trace G G H q (I) such that the transmission problem 

(A/^ - §)$,ef = i^ief in Hlef, $i«f(0,t) = G(t) + G°(t) + (co(t), 

« g - |)$ rig = F rig in Hrig, $ rig (0,t) = G(t), 

dA f {0, t) - <9 CT $ rig (0, t) = H{t) on /, 
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admits a (unique) solution ($| e f, $ r ig) with exponential decay. 

Proof. Let ($® ef , Co) be the solution provided by Lemma 5.2 for the data F = F\ e f and G = 0. 
Let be the unique exponentially decreasing solution of the problem 

(A/7-i)$° ig = F rig in Hrig, $° g (0,t) = 0. 

Let H° be the jump 9 CT $° ig (0, £) -^$^(0, t). If we define the new unknowns $^ ig = $ rig -$r ig 
and $l ef = $| ef — $° ef , the problem we want to solve becomes 

(^-1)^ = in Hlef, $ f l ef (0,t) = G(t) + (C-Co)co(t), 

1)^ = in Hrig, 3^(0,*) =G(t), 
8^(0,1) -d a $l ef (0,t) = H(t)-H°(t) on /. 

Using Proposition 6.3 we can set G = (T ng +T lef il 1 )~ 1 (i/— H ), which ensures the solvability 
of the above problem. □ 

6.1.2. Construction of quasimodes. 

• Terms of order h°. Let us write the "interior" equations: 



lef. 



£()^lef,0 = 70*lef,0 



lef* : A/" lef $ief, = 70$lef,0 

rig : K g $ rig ,o = 7o$ri g ,o • 

The boundary conditions are: 

*lef,o(0,t) + $|ef,o(0,t) = $rig,o(0,t), 
^$|ef, (0,t) = 9 CT $ ri g,o(0,t). 

We get: 

7o = r, *ief,o = 5 , o(s)c (t). 
o 

We now apply Lemma 6.4 with i 7 ^ = 0, F rig = 0, G = 0, H = to get 

G = and £ = 0. 

We deduce: $i e f,o = 0, $ r ig,o = and, since £ = — <7o(0), <?o(0) — 0- At this step, we do not 
have determined g yet. 

• Terms of order h 1 ^ 3 . The interior equations read: 

lef s ■ £ W|ef,i = 70^lef,l + 7l*lef,0 

lef ff : A^ ef $, efil = 7o$ief,i + 7i$ief,o 

rig : Ar o ng $ rig ,i = 7o$ri g ,i + 7i$ri g ,o- 

Using Lemma 5.3, the first equation implies: 

7i = 0, #itf,i(s, t) = gi(s)co(t). 
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The boundary conditions are: 

0l(O)co(t) + $|ef,l(O,t) = $ rig ,l(0,t), 

g' (0)c Q (t) + 5 ff $ief,i(0,t) = 5 CT $ rig! i(0,t). 

The system becomes: 

lef CT : (A^ ef -i)$ lefil = 0, 

ng : (K 6 - ^K,i = 0. 

We apply Lemma 6.4 with F| ef = 0, F rig = 0, G = 0, = -g' (0)c (t) to get: 

C7 = -^(0)(T ri s + T lef n 1 )- 1 c . 
Since G = $ r i g ,i and C = — «7i(0), this determines $i e f,i, Q r \g,i and <7i(0). 
• Terms of order h 2 ^ 3 . The interior equations write: 

lef s : C 2 ^\ef,0 + £o*lef,2 = Tl^lrf,*; 

i+fc=2 

lef CT : A^ ef $lef,2 = ^ 7/$lef,fc 

l+k=2 

rig : K g $ rigi2 = i$ rigi2 , 



with 



£2*lef,0 = -g'o(s)co(t) + — -=S^o(s)5 t 2 (c ). 

Lemma 5.3 and then Lemma 5.4 imply: 

(6-10) - g'o - ^-^ s 9o = l29o- 

Thus, 72 is one of the eigenvalues of the Airy operator and g an associated eigenfunction. In 
particular, this determines the unknown functions of the previous steps. We are led to take: 

#lef,2M) = *irf >2 + 92(s)co(t), With tfj^ 2 = 



and to the system: 



lef CT : (A^ ef -i)$ lef , 2 = 



rig : (K g - ^)<f rig ,2 = 0. 

Using Lemma 6.4, we find 

C7 = -^ 1 (0)(T ri s + T lef n 1 )- 1 Co . 
This determines $ r i g ,2, $i e f,2 and #2(0). The function g\ is still unknown at this step. 
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• Further terms. Let us assume that we can write \&ief,fc = & + gk{ s ) c o{t) for < A; < 
j and that (gk)o<k<j-2 and (\l/j^ f k )o<k<j are determined. Let us also assume that g.,_i(0), 
(7fc)o<fc<j, ($ri g ,fe)o<fc<j-i, ($ief,fc)o<fc<j-i are already known. Finally, we assume that ^(0), 
$ief ,j, ^rigj are known once g,j^\ is determined and that all the functions have an exponential 
decay. 

Let us collect the terms of order h^ +1 ^ 3 . The interior equations write: 

i+i j'+i 

fc=0 fc=0 

j'+i j+i 

\ef a : A^ ef $| ef , i+ i_ fc = 53 7fe$|ef j+i-fc 

fc=0 fc=0 

i+i j+i 
rig : A^ lg $ rigii+ i_ fc = ^ 7fc^rig,i+i-fc, 

fc=0 fc=0 

We examine the first equation and notice that C\ = and 71 = so that ^i e f j does not appear. 
We can write this equation in the form: 

£o - g j *lef,j+l = -A*lef,i-1 - 72*lef,i-l ~ 7j+l*lef,0 

— 53 A^lefJ+l-fc ~~ 5Z Tfe^lef j'+l-fe- 
fe=4 fc=3 

We apply Lemma 5.3 and we obtain an equation in the form: 

1 

where / and ^_i(0) are known. Then, Lemma 5.4 applies and provides a unique value of 7^+1 
such that gfj_! has an exponential decay. From the recursion assumption, we deduce that gj(0), 
®\ef,j, ^rigj are now determined. Lemma 5.3 uniquely determines \&j^ f - +1 such that: 

We can now write the system in the form: 

lef CT : (A^ ef -l)$ lefj+1 = F lef 

where i*| ef , F rig have an exponential decay. The transmission conditions are, cf. (6.7)-(6.8): 

$lef,i+l(0,t) = $ rigJ+ i(0,t) - *| efj+ l(0,t) 

= $rigj+i(0,0 - ^ fj+ i(0,t) - <7 j+ i(0)cb(i) 

and 

9 a $i efJ+1 (o,t) - a CT $ rig> i+i(o,0 = # (t) = -^(o)co(t) + H(t), 
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where H is known. We can apply Lemma 6.4 which determines $ ri g,j+i> ^lefj+i (with an 
exponential decay) and <?j+i(0) once gj is known. 



Quasimodes. The previous construction leads to introduce: 

( J+2 



(6.11a) #W) = { 



lefj 



^ j=0 

where the correction term 

(6.11b) Rj, h (r) = d s ^ ]e{J+2 (0,T)h J / 3 

J+2 



h" 3 



when u < 



when m > 



E 



td t 



7T 



\ J + 2 a 

(0, r) + $ lefj (0, r)) J + ^ ^*W°> r) ^ /3 



lefj 



is added to make ^ satisfy the transmission condition (6.8). Here we have used a smooth 
cutoff function x" e being 1 near 0. By construction, ip^ defined by the identity 

(6.11c) ^ ] (x,y) = X ]ef (u)^ ] (u,-) 

belongs to the domain of Cq u \(K). Using the exponential decays, for all J G N we get the 
existence of h > 0, C(J, h ) > such that for h G (0, ho): 

J+2 

(£ Gui (h) -J2^ h3/3 )^h ] \\ ^ C ^ M h '~ 

j=0 



L+J/3 



6.2. Agmon estimates and consequences. In this last subsection, we prove Theorem 2.6. For 
that purpose, we first state Agmon estimates to show that the first eigenfunctions are essentially 
living in the triangle Tri so that we can compare the problem in the whole guide with the 
triangle (see also Section 4.4 where this idea was explained in the one-dimensional setting). 

Proposition 6.5. Let (\,ij)) be an eigenpair of Cq u \(K) such that |A — ~| < Ch 2 ^ 3 . There exist 
a > 0, ho > and C > such that for all h G (0, ho), we have: 



x>0 



hd x ip\ 2 ) dxdy < C 



Proof. The proof is left to the reader, the main ingredients being the IMS formula and the fact 
that "HecGui is a lower bound of £Gui(^) hi the sense of quadratic forms (see the analysis of 
Propositions 4.3 and 4.4). See also [14, Proposition 6.1] for a more direct method. □ 

• Proof of Theorem 2.6. Let ijj^ be an eigenfunction associated with AG U i,n(^) and assume 
that the ip^ are orthogonal in L 2 (Vt), and thus for the bilinear form B Gu - ljh associated with the 
operator £ Gui (/i). 

We choose e G (0, |) and introduce a smooth cutoff x h at the scale h l ~ £ for positive x 



e-l> 



with \ '. = 1 if x < |> X — if x > 1 



PLANE WAVEGUIDES WITH CORNERS IN THE SMALL ANGLE LIMIT 

and we consider the functions x^n- We denote: 
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<^ o (/0=span( X ViV--,xVX)- 

We have: 

and deduce by the Agmon estimates of Proposition 6.5: 

QG U i,UxVn) = [\GuUh) + 0{h^))\\ X h ^ h n \\ 2 . 

In the same way, we get the "almost"-orthogonality, for n ^ m: 

iW(x%xNt) = ^°°)- 

We deduce, for all v G <£ No (h): 

Qgu;M < (^,N a (h) + 0(h™))\\v\\ 2 . 

We can extend the elements of <£ No (h) by zero so that QGu\,h(v) = Qjri eth (v) for v £ (E No (h) 
where Tri^ is the triangle with vertices (— ny/2, 0), (/i 1_e ,0) and (h l ~ £ ,h l ~ £ + Try/2). A 
dilation reduces us to: 

on the triangle Tri. The lowest eigenvalues of this new operator admits the lower bounds 
| + 2; A (n)/i 2 / 3 — Ch 1 ' 6 ; in particular, we deduce: 

X Gui , No (h) > \ + ^(N )h 2 ^ - Ch 1 ^. 

o 

This provides the separation of the eigenvalues and, joint with Proposition 6.1, this implies 
Theorem 2.6. 
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6.3. Conclusion. 

6.3.1. Eigenfunction asymptotics. With Theorem 2.6, we deduce that the lowest eigenvalues 
of Cq u \{K) are simple as soon as h is small enough. Then, through the spectral theorem, we 
infer that the quasimodes constructed in (6.11) are approximations of the true eigenf unctions 
(see for instance [24]). As a consequence, with the coordinates u,t,r defined in (5.1) and 
(6.2), the n-th normalized eigenfunction admits the following expansion: 

[Ei>o (*n,ief,i(p73,*) + $n,i efj (f,t)) h^ 3 when u < 

(6J2) ^ } ~ £ *W (I r) h>" when u > , 

I i>o 

where the functions ^ n ,ief,j> 3>n,iefj> 3?7i,rigj were constructed in Section 6.1.2 (the subscript 
n emphasizes the dependence on the rank of the zero of the Airy function determined when 
solving Equation (6.10)). 

6.3.2. Remark on the Born-Oppenheimer approximation. At the very beginning of this paper 
we have introduced the operator V.BO,Gu\(h) (see (2.8)) and we have somehow suggested that 
it is an approximation of £Gui(^) in the limit h — > 0. It turns out that we have not used 
7~(-BO,Gu\(h) to investigate the spectrum of the waveguide. In fact, our analysis proves that 
the two term asymptotic expansion of the eigenvalues of V.BO,Tr\(h), "HBO,Gui(^)> Cj r \{h) and 
CGui(h) are the same so that we can a posteriori say that V.BO,Gu\(h) approximates Cc u ](h). 

6.3.3. Back to the physical coordinates. The two-term asymptotics 

*¥(s,t)l 8<0 + h 1 / 3 ^ (a,t)l CT<0 + $? g (<7,T)l ff > ) 

provides us with the leading behavior of the eigenvectors in the scaled half-guide fl It is 
interesting to come back to the physical domain and to interpret this two-term asymptotics 
in the original variables (xi,X2). We have to chain formulas (2.2) giving (x, y), (5.1) giving 
(u, t), (6.2) giving (u, r), and (5.5) giving s and a. We have also to take the relation h = tan 9 
into account. 

Returning to section 6.1.2 and more particularly to (6.10) — and Lemma 5.4, we find that 

< (M) = A((4vrv / 2)" 1/3 S + * A (n)) cos (y). 
Coming back to physical variables (xi, x 2 ) we find that 

<( S) t) = A((|-) 1/3 x 1 + ^ A H)cos(^-^)+^ 2 ) as 6^0. 

As for the term $i := $ l L ef l (7< o + $ r 1 ' s l (7 >o, we find that there exists a profile $i independent 
of 9 such that 



$i(<r,tla<o + r1U>o) = $i(xi,x 2 ) + C>(# 2 ) as 9 -> 0. 



PLANE WAVEGUIDES WITH CORNERS IN THE SMALL ANGLE LIMIT 



37 



Here Xi = x 1 and 

7TX 2 COS 9 



if xi < 0, 



X2 



Ti + X\ sin 9 

X2 cos 9 — x\ sin 9 if x\ > 0. 

This profile $i is exponentially decreasing as x± — > ±oo. It is solution of a transmission 
problem with smooth data for the Laplace operator on the infinite strip M. x (0, ir) with mixed 
Neumann-Dirichlet conditions on the bottom side x 2 = 0, and Dirichlet on x 2 = n. Hence, it 
is piecewise if 2 modulo the addition of a multiple of the singular function V'sing' °f- (1-2). 

The consequence of this is that the coefficient of the singularity f° r a normalized eigen- 
vector of AgJ behaves as 0{9 1 l z ) as 9 ->■ 0. 

6.3.4. X-shaped waveguides. Our results provide without any difficulty the structure of the 
eigenpairs of lowest energy in the small angle limit when the domain is formed by the union 
of two infinite strips of same width ix crossing with an angle 29 (this model appears in the 
physical literature, see [8]). The two non-convex corners of this structure are at the dis- 
tance -P-Q = 0{9^ 1 ). This X-structure can be viewed as a double symmetric V-structure 
and the eigenmodes can be constructed from the V-structure eigenmodes since they interact 
very weakly (their lower scale is 9 1 ^ 3 ). Nevertheless they do interact by an exponentially small 
tunnelling effect which would be interesting to investigate. 
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